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We investigate the electronic structures of the alloyed Bii-^Sbx compounds based on first-principle 
calculations including spin-orbit coupling (SOC), and calculate the surface states of semi- infinite sys- 
tems using maximally localized Wannier function (MLWF). From the calculated results, we analyze 
the topological nature of Bii-^Sbx, and found the followings: (1) pure Bi crystal is topologically 
trivial; (2) topologically non-trivial phase can be realized by reducing the strength of SOC via Sb 
doping; (3) the indirect bulk band gap, which is crucial to realize the true bulk insulating phase, can 
be enhanced by uniaxial pressure along c axis. (4) The calculated surface states can be compared 
with experimental results, which confirms the topological nature; (5) We predict the spin-resolved 
Fermi surfaces and showed the vortex structures, which should be examined by future experiments. 

PACS numbers: 71.15.Dx, 71.18.+y, 73.20.At, 73.61. Le 



I. INTRODUCTION 

In an ordinary insulator, the valence and conduction 
bands are separated by an energy gap, making it elec- 
trically inert. Therefore, the ordinary insulator is not 
sensitive to the change of boundary condition. Recently 
a new class of insulator, namely topological insulator 
(TI), is proposed^*^. TI also has a bulk energy gap, 
which is usually generated by spin-orbit coupling (SOC); 
however it is different from the ordinary insulator in the 
sense that topologically protected gapless states, robust 
against disorder, appear at the edge or surface of a fi- 
nite sample within the bulk energy gap. Thus the TI 
has conducting channels along its edge or surface. The 
quantum spin Hall (QSH) insulator, such as HgTe/CdTe 
quantum well 2,3,7 , is an example of two-dimensional (2D) 
TI. The conducting edge channels of HgTe/CdTe quan- 
tum wells have been theoretically predicted 2 - and experi- 
mentally observed 3 ' 7 . From the theoretical point of view, 
the TI can be distinguished from the ordinary insula- 
tor by the Z2 topological invariantsA£&£, and the exis- 
tence of gapless spin-filtered edge states on the sample 
boundary is guaranteed for TI. The edge states come in 
Kramers's doublets, and time reversal (TR) symmetry 
ensures the crossing of their energy band at time rever- 
sal invariant momenta (TRIM). Since these band cross- 
ings on the edge are protected by TR, they can not be 
removed by any perturbation respecting the TR symme- 
try, such as non-magnetic impurities. It is expected that 
the robust gapless spin- filtered surface (edge) states have 
novel applications in spintronics. 

Besides the 2D QSH insulator, the TI can also exist in 
three dimensional (3D) materia^£&^±±>±^. Similar to 
the edge states in 2D QSH insulator, in 3D TI, topologi- 
cal surface state protected by TR, which can be described 
by odd number of Dirac points , emerges at the surface of 
the finite 3D sample^^^. Compared with the 2D TI, 



the 3D TI and its surfaces can be readily investigated 
by ARPES and STM experiments. The 3D TI also dis- 
plays the remarkable topological magneto-electric effect 6 . 
Therefore searching for realistic 3D TI is now becoming 
an attractive and challenging subject. 

It was first suggested that the semiconducting alloy of 
bismuth and antimony (Bii-^Sb^) is an example of such 
3D TI 8 ' 11 . Based on the tight-binding (TB) model of Liu 
and Allen 15 , Fu and Kane developed a theory to analyze 
the topological nature of the surface state in Bii-^Sb^ al- 
loy^. Experimentally, Hsieh et aS£ observed the surface 
states by high-momentum-resolution angle-resolved pho- 
toemission spectroscopy (ARPES), and demonstrated 
the topological nature of the surface states by counting 
the number of the Fermi surface crossings from the Zone 
center to the boundary. However, clear discrepancies ex- 
ist between the theory and the experiment about the sur- 
face states, although their final conclusions are consistent 
with each other. On the other side, the surface states of 
pure Bi or Sb have been intensely studied experimentally 
and theoreticall y 16 i 17 i 18 i 19 i 2Q i 21 , but there are still fewer 
careful studies of their alloy. Therefore, in this paper, 
we present a systematic study of the surface states of 
Bii-^Sb^ alloy, based on quantitative first principle cal- 
culations. We show that pure Bi is topologically trivial 
because the SOC is too strong. The effective role of Sb 
doping is to reduce the strength of SOC and revert the 
band ordering at L point of the Brillouin Zone(BZ). Fi- 
nally strong topological insulator can be reached by Sb 
doping. By constructing the maximally localized Wan- 
nier function (MLWF) from the ab-initio schemes, we 
calculate the surface states of semi-infinite system, and 
analyze the shape of Fermi surfaces as well as the spin- 
resolved local density of states. These results are com- 
pared with the experiment of Hsieh et aS£. In addition, 
although direct band gap exists in Bii-^Sb^ system, in- 
direct band gap can be only realized for a very narrow 
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FIG. 1: (color online) (a) Along the (111) direction of rhom- 
bohedral A7 structure, there are three possible atomic posi- 
tions^, B,C) for the triangle plane. Black filled dots denote 
A site, red filled up- triangles denote B site, and blue filled 
down- triangles denote C site, (b) The schematic plot of Bi- 
layers projected onto the plane paralleling to the (111) axis. 
After the dimerization of two Bi layers, the Bi-2 moves to the 
dashed site. Ad denotes the magnitude of the dimerization. 
1st and 2nd nearest neighbor (NN) hopping are the inter- layer 
hopping between the different sublattices, while 3rd NN hop- 
ping is the intra-layer hopping within the same sublattice. 



doping range, and the material has long been regarded as 
typical semi- metallic system. In order to make the bulk 
material insulating completely, based on our calculations, 
we predict that an efficient way to enhance the indirect 
band gap is to apply uniaxial pressure along the c-axis. 

The paper is organized as follows. In sec. [Ill we dis- 
cuss the crystal structure and symmetry of Bii-^Sb^. In 
Sec. mUwe study the transition between the topological 
non-trivial and trivial phases, and present a schematic 
phase diagram as a function of SOC. In Sec. IIVI we de- 
velop an accurate method based on MLWF to obtain the 
surface Green Function, and study the topologically non- 
trivial surface states of Bii-^Sb^. In Sec. [V] we provide 
a brief discussion and conclusion. 



II. STRUCTURE AND SYMMETRY 

Bi and Sb have the same rhombohedral A7 crystal 
structure^ with space group R3m. The A7 structure 
can be regarded as a distorted fee NaCl structure. For 
fee NaCl structure, there are two sets of sublattice, say 
Bi-1 for Na sites and Bi-2 for CI sites, which both form 
the fee structure. In such fee structure, Bi-1 and Bi- 
2 sublattices are equivalent and we can shift Bi-1 sub- 
lattice by (1/2,0,0), or (0,1/2,0), or (0,0,1/2) of cubic 
structure to obtain Bi-2 sublattice. Along (111) direc- 
tion of fee structure, Bi triangle layers are stacked with 
the sequence of ABC ABC ABC • • • , where A, B and C 
denote three different atomic positions for triangle plane, 
as shown in FigHJ Without distortions, there exist two 
kinds of space-inversion center. One is located at the 
Bi layer center and each Bi sublattice is space-inverted 
to itself (called type-I inversion), while the other one is 
located at the middle way of two Bi layers and each Bi 



FIG. 2: (color online) (a) BZ of fee structure, (b) 3D BZ of 
rhombohedral A7 structure and its projection onto the [111] 
surface. The A7 BZ can be obtained from the fee BZ by two 
steps: (1) rotating the (111) direction of cubic to be along the 
c-axis; (2) slightly distortion of the BZ along the c-axis. The 
Li point in fee BZ is changed to be T point in the A7 BZ, 
which is now inequivalent to £2,3,4 due to the distortion. T 
(T), X (L), and K in 3D BZ are projected to T, M and K in 
the 2D BZ of [111] surface. 



sublattice is space-inverted to another sublattice (called 
type-II inversion). 

Starting from fee NaCl structure, two steps are re- 
quired to obtain rhombohedral A7 structure. One is the 
stretching along cubic (111) direction, while the second 
is the relative shift of inter- Bi- layer distances along (111) 
direction, or in other words dimerization of two Bi lay- 
ers. The second step breaks the type-I inversion symme- 
try, but preserve the type-II inversion. Therefore, after 
distortions, the inversion symmetry can only transfer the 
atoms from one sublattice to the other, which is very 
important for our following discussion. 

The Brillouin Zone(BZ) of fee structure and rhombo- 
hedral A7 structure are shown in Fig|2l For fee struc- 
ture, there are four equivalent L points, which are lo- 
cated at Li = (7r,7r,7r), L2 = (7r,— 7r,7r), L3=(7r,7r,— 7r) and 
L4=(7r,— 7T,— 7r) of BZ. After two kinds of distortions 
along cubic (111) direction, the fee structure changes to 
the rhombohedral A7 structure, which breaks the equiv- 
alence between L\ (which is denoted as T point in BZ of 
rhombohedral A7 structure) and £2,3,4- 



III. EFFECT OF ALLOYING AND PHASE 
DIAGRAM 

Although pure Bi and Sb have been studied exten- 
sively, the alloyed system is not carefully considered yet. 
In this section, we will present a simple TB model by tak- 
ing the strength of SOC A and dimerization Ad as two key 
parameters to describe the effect of alloying. Based on 
this TB model, we obtained a schematic phase diagram 
for the topological nature of the compounds, which is in- 
structive for us to understand the main physics. Then 
in the next section, surface states of alloyed system will 
be studied from accurate ab-initio calculations based on 
virtual crystal approximation (VCA). 
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A. Effect of Alloying 

In order to take into account the effect of alloying, cer- 
tain kinds of approximation have to be introduced. The 
first step, which is conventionally followed, is to assume 
the uniform distribution and neglect the disorder effects. 
However, this approximation is not sufficient. To further 
simplify our understanding, we emphasize the following 
factors: (1) Sb is located just on top of Bi in the pe- 
riodical table, therefore they have the same number of 
valence electrons and form the same A7 crystal struc- 
ture. (2) Even 30% Sb alloying into Bi will only modify 
the lattice parameters by around 1%—. (3) The atomic 
SOC strength of Sb is weaker than Bi by a factor of 3. 
Therefore, we believe the strongest effect of Sb alloying 
into Bi is to reduce the SOC strength, and we can neglect 
the effect coming from the change of lattice parameters. 
Following this strategy, we construct a simple TB model 
to understand the main physics. 

We consider one s and three p orbitals of each Bi atoms, 
together with the two sublattices and two spin degree of 
freedom, and totally there are 16 orbitals, denoted as 

\4),\4), bL>> \Pi y ), bL>> \pL), \Ply), \Pu)> \ s i),\4), 

|PL>> \P2y)i IpL)» \P2x)i 1/4,)' \P2z)- Here Z axis is taken 

along the (111) direction of fee structure or the (001) 
direction of hexagonal cell and subscript number denotes 
different sublattice. The hopping parameters are defined 
in FiglU with the expression 



Vppa = V ppa (l + a(Ad - Ad )) 
V pp7T = V pp7T (l + (3(Ad - Ad )) 
V ssa = V ssa (l + j(Ad-Ad )) 
V spa = V spa (l + 5(Ad-Ad )) 
for the first nearest-neighbour (NN) hopping, 

v; p „ = v; p m - p(Ad - Ad )) 

K' sa =K sa (l- 7 (Arf-Ad )) 
V: pa = V; pa (l-5(Ad-Ad )) 

for the second NN hopping and 
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for the third NN hopping. Here the parameters V m , V m 
and are taken from the Liu- Allen modeU^, where m is 
ppa, ppir, ssa or spa, respectively. Since the intra-layer 
Bi-Bi distance is larger than inter-layer Bi-Bi distance, 
V m and V m are larger than V^. Linear dependence of 
Ad due to the dimerization is assumed and Ado is for 
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the experiment structure of pure Bi. When Ad = 0, 
this system has no dimerization, leading to V m = V m . 
Then a set of linear equations is obtained, which can be 
used to determine the value of the parameters a, /?, 7 
and S. Besides the hopping terms, the atomic SOC with 
parameter A is also taken into account. Therefore, the 
final Hamiltonian is given as the function of two variables 
Ad and A, with the form 



(4) 



where B.\\—B.22- l #12=^21 are 8 x 8 matrice. Hu in- 
cludes intra-sublattice hopping and on-site SOC inter- 
action, while H12 represents inter-layer hopping. Since 
the Hamiltonian has the type II inversion symmetry at 
T and L points, a unitary transformation is applied here 
to rewrite the Hamiltonian in the new basis with un- 
ambiguous parity, the odd parity basis, 775 (l^i) — l 5 2))> 

75(Kx) + b2x», + + and 

the even parity basis, ^-(|sf> + |^», ^(K) " K», 
7f(K„} - \P%)), 7f(bL) " \P5z))-_ With the new basis, 
the Hamiltonian is changed to be iJ, 



'Hu 








#22 . 



(5) 



which is block diagonal, because the odd and even par- 
ity states will not mix in a system with space-inversion 
symmetry. 



B. Phase Diagram 

The topological nature of the system can be deter- 
mined from the parity of the occupied bands at TRIM^. 
The band gap for Bii-^Sb^ is near T and L points, there- 
fore here we focus on one T and three L points in BZ. 
The parity of the occupied bands for T and L points can 
be easily obtained since the eigen states of Hu have odd 
parity while those of H22 have even parity. In Fig [3l the 
energy levels of six p bands for both L and T points are 
plotted as a function of Ad, where two different values 
of A are chosen, A = 1.5 corresponding to the value of Bi 
and A = 0.5 corresponding to the value of Sb. Three low- 
est levels of total six p bands should be occupied, namely 
the conduction band and valence band are the third low- 
est band and fourth lowest band respectively, which have 
opposite parities. When increasing the dimerization pa- 
rameter Ad, at T point the band gap increases rapidly 
and there is no leveling crossing between the conduction 
band and valence band, while at L point, the band gap 
is quite small and the sequence of conduction band and 
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FIG. 3: (color online) Energy bands are plotted as a function 
of the dimerization parameter (Ac/) with SOC parameter A 
taken as l.beV and 0.5eV. (a) and (b) are for T point, while 
(c) and (d) for L point. The blue lines denote the states with 
parity +1 while the red ones are the states with parity —1. 
The dashed line represents Ad — Ado • 
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FIG. 4: (color online) The phase diagram of the system with 
two variables: the dimerization parameter Ad and SOC pa- 
rameter A. In the magenta region, the parity of L and T 
points are +1 and -1, respectively, in other words, the sys- 
tem is in the topological non-trivial phase with Z2 invariants 
(1;111). In the green region, however, the parity of L and T 
points are all -1, i.e, it is topological trivial with Z2 invariants 
(0;000). Pure Bi locates in the trivial region and Sb locates 
at topological non-trivial region. The experimental Ad of Bi 
and Sb are almost the same^. 



IV. AB-INITIO CALCULATIONS AND SURFACE 
STATES 



valence band can even change with A = 0.5eV. When 
Ad = Ado, as indicated by the dashed line in Fig 02 
the occupied valence bands for Bi (A = l.beV) and Sb 
(A = 0.5eV) have different parities at L point but the 
same at T point, therefore we conclude the topological 
natures of Bi and Sb are different. This result is the same 
as that of Fu and Kane^, in which they claim that Z2 
invariants (z^o; ^1.^2^3) are (0;000) for Bi, which is topo- 
logically trivial, but (1;111) for Sb, which corresponds to 
topological non-trivial phase. 

By determining the gap closing line of the conduction 
and valence band at L points, we can obtain the phase 
diagram of the alloy Bii-^Sb^ as a function of SOC pa- 
rameter A and dimerization parameter Ad, as shown in 
Fig. [H Since the topological nature of the system can 
only be changed by closing bulk gap, in the magenta re- 
gion, the system should have the same non-trivial topo- 
logical behavior to Sb, while in the green region, it should 
be trivial insulator, which is same to Bi. It can be also 
understood from the parity analysis that the Bi is topo- 
logically trivial because both L and T points have re- 
verted bands ordering, namely parity are -1 for both L 
and T. However, by reducing the SOC strength (i.e. Sb 
doping), the parity of L is recovered to be +1 while T 
remains to be -1, therefore topological non-trivial system 
is realized. In the next section, accurate simulation from 
ab-initio calculations will be presented. 



A. ab-initio Method and Surface Green Function 

The ab-initio calculation is carried out by our BSTATE 
(Beijing Simulation Tool for Atom Technology) 22 code 
with plane wave pseudo-potential method. The gener- 
alized gradient approximation(GGA) of PBE-type 2 ^ is 
used for the exchange-correlation potential. Especially 
in Bi's pseudo-potential, we take into account 5d 10 elec- 
trons as the valence band by ultra-soft pseudo-potential 
scheme. The k-mesh is taken as 12 x 12 x 12 and the 
cut-off energy is 340. OeV for the self-consistent calcula- 
tion. For pure Bi, the optimized lattice parameters are 
a = 4.669l,c = 12.1506A, and d = 0.234lA, which are 
in good agreement with previous calculations. 

To take into account the effect of alloying, virtual crys- 
tal approximation (VCA) is necessary. There are sev- 
eral ways to do this, particularly in the pseudo-potential 
approach, the simple linear combination of Sb and Bi 
pseudo-potentials can be used, and the corresponding 
pseudo-potential of alloyed virtual atom is regenerated 
by solving the atomic problem again^ 2 -. Unfortunately, 
this procedure is not accurate enough for our purpose 
here. As already suggested by previous studies 22 , such 
VCA procedure can be used for those states far away 
from the Fermi level, however for those states very close 
to the Fermi level, the error bar is big. The system we 
study here (Bii-^Sb^) is very sensitive to the p orbitals, 
such VCA pseudo-potential can not give sufficient accu- 
racy. In order to have an accurate VCA scheme, we need 
to consider the particularity of our system. As we already 
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explained in the last section, the alloyed Bii-^Sb^ have 
the same crystal structure and almost the same structure 
parameters. The main effect of Sb alloying is to tune the 
SOC strength, A. Therefore we may have a simple yet 
accurate VCA scheme. Here we take Bi's parameters for 
simplicity and tune A in Bi's pseudo-potential to simu- 
late the doping parameter x of the alloyed Bii-^Sb^. In 
such a way, since we do not need to the solve the atomic 
problem again, the pseudo-potential is accurate enough. 

We are interested in the surface states of the semi- 
infinite system, a method based on maximally local- 
ized Wannier function (MLWF) is developed to calcu- 
late the surface states of semi-infinite system. The ab- 
initio MLW F 24 i 25 method can be regarded as an exact 
TB method with its parameters calculated from ab-initio 
self-consistent electronic structure calculations. First, 
the semi-infinite Bii-^Sb^ system can be divided into two 
parts: the bulk part and the surface part. The bulk part 
Hamiltonian is constructed with the MLWFs from bulk 
Bii-^Sb^ ab-initio calculations, while the surface part 
Hamiltonian is constructed with MLWFs from Bii-^Sb^ 
film with slab calculations. With these MLWFs' hop- 
ping parameters, iterative metho d 26 ! 27 is adopted to solve 
the surface Green function of the semi-infinite system 
G^ Za (k||,e + irj), where n and I denote the super-cell 
along z direction and the atomic bilayer plane within one 
super-cell respectively, a gives the orbital index in one 
atomic bilayer and kn is a good quantum number in semi- 
infinite system. 

The charge density of states (DOS) and spin DOS^ are 
related to the surface Green function with the expression 

<(k|| , e) = - llm G l n a n la (k„ , 6 + «7) (6) 

a 

and 

S l nt<J (k,| , e) = - Ilm £ Gnn W (ky , e + ir,)0^ la 

* a,0 , (7) 

Ola = (iP\s°\la) 

respectively, where s a is the spm(s Xj y iZ ) operator. When 
n = 0, 1 = 0, iVj(k||, e) and S^ a (k\\ , e) give the local DOS 
and local spin DOS at the surface and in the following 
we use iV(k||, e) and ^(kn, e) for short. 



B. Surface State and Fermi surface 

The local DOS N(\t\\,E) at the surface is plotted for 
two different SOC parameters, A=1.28eV and A=l.leV 
for the topologically trivial and non-trivial situations re- 
spectively. In Fig. [5]with A = l.leV, developing from the 
f point, there exists two surface bands connected to the 
M point, which are denoted as Ei and £2 respectively. 
At M point, £2 band returns to valence band while £1 
band merges into the conduction band. Therefore, those 



0.4 
. 3 
§b . 2 

~0 . 1 

8 

UJ 

-0 . 1 

-0.2 

FIG. 5: (color online) Left panel: The surface local DOS for 
A = l.leV. The broad red regions denote the continuous 
bulk bands with a small gap of about lOmeV at M point. 
Two surface bands indicated Si and £2 disperse within the 
bulk gap. The dashed line indicates the Fermi energy, which 
intersects five times with the two surface bands from F to 
M. Right panel: The region framed by the black rectangle in 
the left panel is zoomed in. One surface state (Si) goes up 
to merge into the conduction band while the other one (£2) 
goes back to the valence band. 
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FIG. 6: (color online) Left panel: The surface local DOS for 
A = 1.28eV. The two surface states £1 and £2 disperse in 
a different way in the present case. At the M point, both 
£1 and £2 are connect to the valence band. Therefore the 
Fermi energy intersects four times with the surface bands be- 
tween r and M. Right panel: The region framed by the black 
rectangle in the left panel is zoomed in. 



surface states cross the Fermi energy five times in to- 
tal (odd number), which indicates the topologically non- 
trivial nature of this phase. On the contrary, in Fig. [6] 
(for A=1.28eV), both £1 and £2 bands return to the 
valence band at M point, and they cross the Fermi level 
four times (even number). This indicates that the system 
is topological trivial. 

The shape of the Fermi surface for the two different 
phases is plotted in Fig. [71 and Fig. El There are one f 
point and three M points in the surface BZ, and they 
are all TRIM. The main difference between Fig. [71 and 
Fig. [8] is around the M point. For both cases, the f 
point is enclosed by one Fermi arc, however, the M point 
is different: it is enclosed by one Fermi arc at A = 1.28eV, 
and it is not for A = l.leV. Therefore by counting the 
total number of TRIM enclosed by Fermi surface in the 
BZ, it is even number for A = 1.2SeV (one f plus three 
M), and odd number for A = l.leV (only one f point). 

Here we compare our results with that from TB analy- 
sis^ and that from experiment 10 . As shown in Fig [9] (a), 
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FIG. 7: (color online) Left panel: The Fermi surface plot for 
A = l.leV. The black hexagonal region is the 2D BZ of [111] 
surface for A7 structure. F is enclosed by a hexagonal electron 
pocket. There are other six hole pockets and six electron 
pockets surrounding. Right panel: The region framed by red 
rectangle in the left panel is zoomed in. We can clearly see 
that the outest six small electron pocket don't enclose M. 
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FIG. 8: (color online) Left panel: The Fermi surface plot for 
A = 1.28eV, which is similar to Fig. [7] except that the outest 
six electron pockets enclose the M point. Right panel: The 
region framed by red rectangle in the upper panel is zoomed 



we find five crossing points between the surface bands 
and the Fermi energy along the line from f to M, this is 
the same to those observed in the experiment of Hsieh et 
ulili (Fig [9}, however in a simple TB modelii, the num- 
ber of crossing is three (Fig [9] (b)). A small discrepancy 
is found between our ab-initio calculation and the exper- 
iment of Hsieh et al near M point, as shown in Fig [9] (a) 
and (d). In the experiment of Hsieh et aS£, a third sur- 
face band £3 appears near M point and be degenerate 
with Ei band at M point, however in our calculation, 
there is no such band and £1 band will go up and merge 
with the conduction band. This discrepancy may come 
from additional trivial surface states, as suggested in Fig 
[9] (c). This discrepancy remains to be justified by future 
studies. 

In addition to the energy resolution, we are also able to 
calculate the spin-resolved surface states. As an exam- 
ple, we carry out this calculation for the semi-infinite 
Bii-^Sb^ system's blow surface, and show the spin- 
resolved surface state in Fig. [lOl The spin orientation 
of the surface states at the Fermi level is plotted for 




(a) 



M 
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FIG. 9: Schematic picture for the comparison of the surface 
bands obtained from (a) our ab-initio calculation, (b) TB 
model (from the work of Teo et clIf^) and (d) ARPES ex- 
periment results. In ARPES experiment, an additional £3 
surface band (dotted line in (d)) becomes degenerate with £2 
band at M point. This additional band may come from the 
hybridization between the topological surface states and the 
other trivial surface states, as suggested by the red dotted 
line in (c). 
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FIG. 10: (color online) The spin resolved Fermi surface for the 
semi- infinite Bii-^Sb^'s below surface, whose normal is along 
-z direction. The arrow in (a) indicates (S x ,S y ); different 
colors in (a) and (b) represent S z along different directions. 
The red color means that the S z is along the +z direction; the 
blue color means that the S z is along the -z direction. The 
three pieces of Fermi surface F\ , F2 and F3 are marked in Fig 
[3 A is taken as l.leV here. 
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FIG. 11: (color online) Indirect energy band gap. We de- 
fine the indirect energy gap as the difference between con- 
duction band bottom(CBB) and the valence band top(VBT). 
The blue line and red line present the indirect energy gap's 
functions of c/a for A=1.28eV and A=l.leV. Hi and H2 points 
are not high symmetry points, and locate in the mirror plane 
in BZ. Hi is near to T point and H2 is near to L point. In the 
inset, we show the schematic figures indicating the different 
energy band structures. 

three regions of 2D BZ with A = l.leV. Clearly, ver- 
tex structure is found for the electron pocket around T 
point (-Fi), which confirms the topological nature of sur- 
face state. Because below surface's normal is along the 
-z direction, we can confirm that the chirality of the ver- 
tex structure around T is left-handed and agrees with the 
recent spin-resolved ARPES experimen t 2829 . 

C. Indirect Band Gap 

In the above discussions, we call Bii-^Sb^ bulk as "in- 
sulator" because there exists a direct band gap between 
the conduction and the valence bands. Unfortunately, 
Bii-^Sb^ is actually a semi- metal (not true insulator) 
for most of the doping range x, namely there exists finite 
overlap between the conduction band bottom (CBB) and 
the valence band top (VBT). If we define the true gap E g 
(indirect gap) as the energy difference between the CBB 
and the VBT, E g is negative for most of the x, and it is 
positive only for 0.07< x <0.22 n . On the other hand, 
to identify the TI nature, except the parity arguments as 
discussed above, it is crucially important to have a full 
bulk gap throughout the BZ. Therefore, a serious ques- 
tion for Bii-^Sb^ is "can we make the indirect gap E g as 
positive as possible?" , or in other words, "can we widen 
the range of doping x where system is truly insulating?" . 
Here we will show that applying the uniaxial pressure is 



an efficient way to open up the indirect band gap E g . 

Fig. 11 shows the calculated indirect band gap E g as 
function of c/a ratio with fixed volume. The c/a ratio 
can be tuned either by c-axis pressure or by forming thin- 
film matched to substrate with different lattice parame- 
ters. For both sides of the topological phases (A=1.28eV 
or l.leV), a broad positive E g region can be obtained by 
reducing c/a ratio slightly (around 3% reduction from 
its experimental value c/a=2.6). For both A=1.28eV 
and A=l.leV, the CBB is located at L point, however 
the VBT is located at T point for A=1.28eV, and at Hi 
point (around T point in the mirror plane) for A=l.leV. 
Despite of the different positions of VBT, the effect of 
c-axis pressure is always to raise the energy levels round 
L point, and lower the levels around the T point. For 
A=l.leV, energy level at Hi point goes lower, and en- 
ergy level at H2 point near the L point in the mirror 
plane goes upper. Therefore the positive indirect band 
gap is realized as schematically illustrated in the insets 
of Fig.ll. 



V. CONCLUSIONS 

As a summary, we develop a method to study the al- 
loyed Bii-^Sb^ system and present a phase diagram to 
describe the topological nature of the system. We show 
that Bi is topologically trivial because the SOC is too 
strong. By alloying with Sb, the effective SOC strength 
is reduced and the topologically non-trivial phase is real- 
ized. By accurate ab-initio calculations and MLWF, we 
calculate the surface states of semi-infinite system. The 
results are compared with recent experiments. We pre- 
dict the spin-resolved Fermi surface which can be tested 
by spin resolved ARPES. Finally, we suggest an efficient 
way to tune the indirect band gap by uniaxial pressure, 
such that true bulk insulating state can be realized for a 
broad doping range. 
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